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Abstract 

We have used numerical micromagnetics for the calculation of the magnetic (small-angle) neu- 
tron scattering cross section of nano composites. The novel aspect of our approach consists in the 
possibility to study the applied-field dependence of the individual contributions to the total mag- 
netic scattering. Such a micromagnetic tool ideally complements neutron experiments in which 
one generally measures only a weighted sum of the Fourier components of the magnetization. The 
procedure furnishes unique and fundamental information regarding the magnetic microstructure 
and corresponding magnetic scattering from nanomagnets. In particular, our simulation results 
explain the recent observation of dipolar correlations in two-phase nanocomposites and provide an 
answer to the question of the explicit dependence of the magnetization Fourier coefficients on the 
scattering vector. 
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Introduction. — Magnetic neutron scattering and, in particular, magnetic small-angle neu- 
tron scattering (SANS) is a very powerful technique for the investigation of spin structures 
in magnetic materials on a length scale between ~ 1 nm and a few hundred of nm (for 
recent reviews see, e.g., Refs. [l-3]). The unique feature of SANS is the possibility to study 
magnetic structure in the bulk of materials, in contrast to various spectroscopic techniques, 
which mostly provide information about the magnetization state at or near the sample sur- 
face. Magnetic SANS has been employed, for instance, to study the range of magnetic 
correlations in nanocrystalline 3d transition metals j^, 5h the vortex lattice of type-II su- 



perconductors 



6|], magnetization dynamics in ferroffuids 



permanent magnets 



7], magnetic domains in Nd 2 Fei 4 B 



8j, nanocrystalline Tb with random paramagnetic suceptibility [9j], the 



spin- helix chirality in FeCoSi single crystals 10J, so-ca 
field-induced mag netization in multiferroic HoMnC>3 



led skyrmions in MnSi 



electric- 



12], the spin structure of core-shell 



13] , or the impact of heterogeneities on the magnetostriction of FeGa alloys 



nanoparticles 

Q- 

The quantity of interest in a magnetic SANS experiment is the elastic magnetic dif- 
ferential scattering cross section dH/dQ, which is usually recorded on a two-dimensional 
position-sensitive detector. Basic scattering theory prescribes that dH/dQ can be expressed 
in terms of the Fourier coefficient M = M(q) of the magnetization vector M(x), more 
specifically, dH/dVL is a weighted sum of the products of Cartesian components of M. For 
bulk ferromagnets, M depends in a complicated manner on the momentum-transfer vector, 
applied magnetic field, and on the magnetic interaction parameters (exchange interaction, 
magnetic anisotropy, dipolar interaction), and only for special cases, e.g., in the approach- 
to-saturation regime, one can obtain approximate closed-form expressions for M [2]. 

The fact that the experimental SANS pattern dH/dQ is composed of several individ- 
ual contributions often hampers the straightforward interpretation of recorded SANS data. 
While, in principle, some Fourier coefficients are accessible by the experiment, e.g., through 
the application of a saturating magnetic field or by exploiting the neutron-polarization de- 



15]). it is often 



gree of freedom via so-called SANSPOL or POLARIS methods (e.g., jsl, Q, 
difficult to unambiguously determine a particular scattering contribution without "contam- 
ination" by unwanted Fourier components. For instance, when the applied field is not large 
enough to completely saturate the sample, then the scattering along the field direction does 



not represent the pure nuclear SANS, but contains also the magnetic SANS due to the 



misaligned spins 



In this Letter, we report the results of full-scale three-dimensional micromagnetic sim- 
ulations of the magnetic SANS cross section of magnetic nanocomposites. Both numerical 
micromagnetics [17| and magnetic neutron scattering are well developed and established 
methods which are widely employed for studying magnetism in solid-state physics. As we 
will show in the following, it is their combination which provides new insights into the fun- 
damentals of magnetic SANS and, thus, into the magnetic microstructure of nanomagnets. 
In particular, the decisive advantage of this approach resides in the possibilitiy to study 
the contributions of the individual Fourier components of the magnetization to dTj/dVl — 
rather than their combination — and relate them to the underlying magnetic microstructure. 
This sheds new light on the ongoing discussion regarding the explicit momentum-transfer 
dependence of dT^/dfl [18]. The micromagnetic computations have been adapted to the 
microstructure of the iron-based two-phase alloy NANOPERM for which experimental data 
exist [19 ]. 



Details of the micromagnetic algorithm. — In our micromagnetic model we have taken 
into account the four standard contributions to the total magnetic energy: external field, 
(uniaxial) magnetic anisotropy, exchange and dipolar interaction energies. The two-phase 
nanocomposite microstructure, consisting of magnetically "hard" iron-based particles em- 
bedded in a magnetically "soft" amorphous matrix, was generated by employing an algo- 
rithm described in Ref. [20]. The simulation volume (= sample volume) is a rectangular box 
of size 125x380x380 nm 3 , which was discretized into N = 10 5 mesh elements. The average 



size of a "hard" inclusion (nanocrystal) is D = 10 nm (as in NANOPERM 19j), whereas 
the mesh size used to discretize the "soft" phase is two times smaller. This discretization 
scheme then limits the accessible range of momentum transfers (via the sampling theorem) to 
q < g max = 1 nm -1 . The volume fraction of the nanocrystallites is Xq = 40 %, correspond- 
ing to about 8000 nanocrystals in the simulation volume. Materials parameters for hard 
("/i") and soft ("s") phases are: magnetizations = 1750 kA/m and M s = 550 kA/m, 
anisotropy constants Kh = 4.6 x 10 4 J/m 3 and K s = 1.0 x 10 2 J/m 3 . As a value for the 
exchange-stiffness constant we used A = 0.5 x 10~ n J/m for interactions both within the 
soft phase and between the hard and soft phases. Equilibrium magnetization state of the 
system was found, as usual, by minimizing the total magnetic energy (for more details on 



our micromagnetic methodology see Ref. 20]). The computed SANS cross sections shown 



below represent averages over typically 8—16 independent random configurations of the hard 
crystallites. 

Magnetic SANS cross section. — For the most commonly used scattering geometry in a 
magnetic SANS experiment, where the applied magnetic field H || e z is perpendicular to 
the wave vector ko || e x of the incident neutrons, the elastic magnetic SANS cross section 
dTj/dVt for unpolarized neutrons can be expressed as 

f^(q) = (\M X \ 2 + \M y \ 2 cos 2 9 + |M| 2 sin 2 9 

— (M y M* + M*M Z ) sin 9 cos 9 



V is the scattering volume, bu = 2.699 x 1CT 15 m//ie (f^B- Bohr magneton), c* is a quantity 
complex-conjugated to c, 9 denotes the angle between the scattering vector q and H, and 
^0,y,z)(q) are the Fourier transforms of the magnetization components M( Xi2/j2 )(x). Note 
that in the small-angle limit and for this particular geometry q = q (0, sin 9, cos 9). Since the 
focus here is on magnetic spin- misalignment scattering, we have ignored the nuclear SANS. 

Results and discussion. — Figure [T] displays projections of the functions jM^ 2 , |My| 2 , 
|M Z | 2 , and of the cross term CT := —(M y M* + M*M z ) into the plane of the two-dimensional 
detector at selected external-field values. Figure [2] shows the field dependence of the mag- 
netic SANS cross section dY^/dVL (computed by means of the above expression for dTj/dVt) 
and of the so-called difference cross section, where dT, / dQ at complete saturation (upper row 
left image in Fig. [2]) has been subtracted from the cross section at the respective field. It can 
be seen in Fig. □ that both \M X \ 2 and \M Z \ 2 are isotropic (i.e., 9 independent) over the whole 
field and g-range. By contrast, at the smallest q and largest fields, the Fourier coefficient 
| My | 2 reveals a pronounced angular anisotropy with maxima roughly along the diagonals 
of the detector (the so-called "clover- leaf" anisotropy), whereas at the smaller fields, the 
anisotropy of \M y \ 2 is rather of the cos 2 #-type (i.e., elongated parallel to H). At saturation 
(fj>oH = 5 T), both | Ma- 1 2 and |My| 2 are relatively small and the main contribution to dYi/dVt 
is due the term |M 2 | 2 , which originates from nanoscale jumps of the magnetization at phase 
boundaries. On decreasing the field, the transversal components increase in magnitude as 
long-range spin misalignment develops. The CT oscillates in sign between quadrants on the 
detector; it is positive for 0° < 9 < 90°, negative for 90° < 9 < 180°, and so on. When the 
CT is multiplied by sin 9 cos 9, the corresponding contribution to the total dT^/dfl becomes 



positive-definite for all angles 9. Therefore, and contrary to the common assumption that 
the CT averages to zero for statistically isotropic polycrystalline microstructures, the CT 
appears to be of special relevance in nanocomposite magnets. 

The finding that \M X \ 2 and \M Z \ 2 are isotropic and that \M y \ 2 = \M y \ 2 (9) provides a 
straightforward explanation for the experimental observation of the clover-leaf anisotropy 



in the SANS data of the alloy NANOPERM 19]. Our simulation results for the difference 



cross section oc (|M r | 2 +|M y | 2 cos 2 ^-l-CTsin^cos^) (see Fig. [2]) agree qualitatively well with 
the experimental data [20J. Note also that clover-leaf-type anisotropics in dTj/dVt have been 



reported for a number of other materials, including precipitates in steels 
gadolinium 



16| . nano crystalline 



2l| . and nanoporous iron [22]. The maxima in \M y \ 2 depend on q and H, and 
on the magnetic interaction parameters and may appear at angles 9 significantly smaller 
than 45°, e.g., at 6 = ±30° (compare Fig. [3}. 

As qualitatively discussed in Ref. Jill, the appearance of the clover- leaf anisotropy in 
dYi/dTt is related to the particular 9 dependence of M y , which is imparted by virtue of the 



magnetodipolar interaction 23j. In fact, up to now, the physical origin for the existence 
of the clover-leaf in the magnetic SANS cross section was merely discussed in relation to 
the jump AM in the magnetization magnitude at the interface between the Fe particle 



and the amorphous magnetic matrix (AM ^ 1200 kA/m for NANOPERM [19j]). This 
jump in magnetization gives rise to an inhomogeneous magnetodipolar field which decorates 
each nanoparticle and which causes nanoscale spin deviations within the matrix in the 
vicinity of each nanoparticle. As an illustration, Fig. H] displays the real-space magnetization 
distribution around two nanoparticles. The symmetry of the spin structure replicates the 
symmetry of the CT (compare to Fig. [T|). In the presence of an applied magnetic field 
the stray-field and associated magnetization configuration around each nanoparticle "look" 
similar (on the average), thus giving rise to dipolar correlations which add up to a positive- 
definite CT contribution to the magnetic dH/dQ. 

Next, we demonstrate that magnetodipolar correlations and the corresponding contri- 
bution of the CT to the magnetic SANS cross section are of relevance for practically all 
bulk magnetic materials which exhibit spatial variations in the magnetic parameters. In 
particular, not only variations in the magnetization magnitude (and possibly exchange cou- 
pling), but also variations in direction and/or magnitude of magnetic anisotropy K (ran- 
dom anisotropy) may give rise to corresponding dipolar correlations. In order to study the 



impact of such variations in K (which are, by construction, naturally included into our 
micromagnetic algorithm), we have computed the spin distribution for the situation that 
M h = M s = M (i.e., AM = 0) but for different values of M. The results for \M y \ 2 are 
summarized in Fig. Note that \M X \ 2 and \M Z \ 2 are both isotropic in this case (data not 
shown). 

Figure |5] reveals that a clover- leaf- type pattern in \M y \ 2 develops with increasing magne- 
tization value M, i.e., with increasing strength of the magnetodipolar interaction. As jumps 
in M at phase boundaries are excluded here as possible sources for perturbations in the 
spin structure, it is straightforward to conclude that nanoscale fluctuations in K give rise 
to inhomogeneous magnetization states (with V ■ M 7^ 0), which decorate each nanoparticle 
and which look similar to the structure shown in Fig. HJ This observation strongly sug- 
gests that the origin of the clover- leaf pattern in dH / dfl of nanomagnets is not only related 
to variations in magnetization magnitude but also due to variations in the magnitude and 
direction of the magnetic anisotropy field. 

Summary and conclusions. — Using a recently developed micromagnetic simulation method- 
ology we have computed the magnetic small-angle neutron scattering (SANS) cross section 
dYi/dVL of a two-phase nanocomposite. This approach allows one to study the applied-field 
dependence of the individual scattering contributions to dT,/dQ, in contrast to experiment, 
where generally only a weighted combination of the magnetization Fourier coefficients is 
measured. It is this particular circumstance, in conjunction with the flexibility of our mi- 
cromagnetic package in terms of microstructure variation (particle size and distribution, 
materials parameters, texture, etc.), which makes us believe that the approach of combining 
full-scale three-dimensional micromagnetic simulations with experimental magnetic-field- 
dependent SANS data will provide fundamental insights into the magnetic SANS of a wide 
range of magnetic materials. As we have demonstrated for the example of the iron-based 
two-phase alloy NANOPERM, we were able to explain on a deeper level the physical origin 
of the recently observed clover-leaf angular anisotropy in the magnetic SANS cross section. 
As a quite general result, our micromagnetic simulations suggest that magnetodipolar cor- 
relations (and the associated clover- leaf-shaped pattern in dH/dQ) are of relevance for all 
bulk magnets with inhomogeneous magnetic interaction parameters. 
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FIG. 1. (color online). Results of the micromagnetic simulations for the Fourier coefficients of 
the magnetization. The images represent projections of the respective functions into the plane 
of the detector (i.e., q x = 0). The external magnetic field H || e 2 is applied horizontally in the 
plane of the detector. Values of H decrease from top row (5 T) to bottom row (30 mT) (see 
insets). From left column to right column: \M X \ 2 , \M y \ 2 , \M Z \ 2 , and CT = —(M y M* + M*M Z ). 
Materials parameters of NANOPERM were used (see text). Pixels in the corners of the images 
have q = 0.8 nm" 1 , Logarithmic color scale is used. In the first three columns from left, red 
color corresponds to "high" and blue color to "low intensity"; in the fourth column, blue color 
corresponds to negative and red color to positive values of the CT. 



5000 mT 300 mT 100 mT 30 mT 




FIG. 2. (color online). Applied-field dependence of the total magnetic SANS cross section dT,/dQ 
(upper row) and of the difference cross section (lower row) (where the SANS cross section at 
a saturating field of hqH = 5 T has been subtracted from the dS/dQ at the respective field). 
Materials parameters and all other settings are as in Fig. [TJ 
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FIG. 4. (color online). Two-dimensional real-space image of the computed spin distribution around 
two nanoparticles (violet circles). Red arrows: magnetization component Mj_ perpendicular to the 
applied field (H is horizontal in the plane); thickness of arrows is proportional to the magnitude 
of Mj_. Blue lines: dipolar field distribution. 



250 kA/m 500 kA/m 750 kA/m 1030 kA/m 




FIG. 5. (color online). Fourier coefficient \M y \ 2 at n H = 0.3 T and for M h = M s = M (i.e., 
AM = 0). M increases from left to right (see insets). Kh = 4.6 x 10 4 J/m 3 , K s = 1.0 x 10 2 J/m 3 
and random variations in easy-axis directions from particle to particle. All other settings are as in 

Fig.m 
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